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Abstract We are living in a dusty universe: dust is ubiquitously seen in a wide variety of 
astrophysical environments, ranging from circumstellar envelopes around cool 
red giants to supernova ejecta, from diffuse and dense interstellar clouds and 
star-forming regions to debris disks around main-sequence stars, from comets to 
interplanetary space to distant galaxies and quasars. 

These grains, spanning a wide range of sizes from a few angstroms to a few 
micrometers, play a vital role in the evolution of galaxies as an absorber, scat¬ 
tered and emitter of electromagnetic radiation, as a driver for the mass loss of 
evolved stars, as an essential participant in the star and planet formation process, 
as an efficient catalyst for the formation of H 2 and other simple molecules as 
well as complex organic molecules which may lead to the origins of life, as a 
photoelectric heating agent for the interstellar gas, and as an agent shaping the 
spectral appearance of dusty systems such as protostars, young stellar objects, 
evolved stars and galaxies. 

In this review I focus on the dust grains in the space between stars (interstellar 
dust), with particular emphasis on the extinction (absorption plus scattering) and 
emission properties of cold submicron-sized “classical” grains which, in thermal 
equilibrium with the ambient interstellar radiation field, obtain a steady-state 
temperature of ~ 15-25 K, warm nano-sized (or smaller) “ultrasmall” grains 
which are, upon absorption of an energetic photon, transiently heated to temper¬ 
atures as high as a few hundred to over 1000 K, and the possible existence of a 
population of very cold (< 10 K) dust. Whether dust grains can really get down to 
“temperature” less than the 2.7 K cosmic microwave background radiation tem¬ 
perature will also be discussed. The robustness of the silicate-graphite-PAHs 
interstellar dust model is demonstrated by showing that the infrared emission 
predicted from this model closely matches that observed for the Milky Way, the 
Small Magellanic Cloud, and the ringed Sb galaxy NGC 7331. 
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1. Introduction 

The space between the stars (interstellar space) of the Milky Way Galaxy 
and other galaxies is filled with gaseous ions, atoms, molecules (interstellar 
gas) and tiny dust grains (interstellar dust). The first direct evidence which 
pointed to the existence of interstellar gas came from the ground-based de¬ 
tection of Na and Ca+ optical absorption lines (Hartmann 1904).^ This did 
not gain wide acceptance until Struve (1929) showed that the strength of the 
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Ca+ K-line was correlated with the distance of the star. The true interstel¬ 
lar nature of this gas was further supported by the detection of the first in¬ 
terstellar molecules CH, CH+ and CN (Swings & Rosenfeld 1937, McKellar 
1940, Douglas & Herzberg 1941). The presence of dark, obscuring matter 
in the Milky Way Galaxy was also recognized at the beginning of the 20th 
century (e.g. see Barnard 1919). Trumpler (1930) first convincingly showed 
that this “obscuring matter” which dims and reddens starlight consists of small 
solid dust grains. The dust and gas are generally well mixed in the inter¬ 
stellar medium (ISM), as demonstrated observationally by the reasonably uni¬ 
form correlation in the diffuse ISM between the hydrogen column density A^h 
and the dust extinction color excess or reddening E{B — V) = As — Ay- 
Nii/E{B — F) 5.8 X 10^^ mag“^ cm“^ (Bohlin et al. 1978), where As 
and Ay are the extinction at the B (A = 4400 A) and V (A = 5500 A) wave¬ 
length bands. From this correlation one can estimate the gas-to-dust ratio to 
be ~210 in the diffuse ISM.^ Despite its tiny contribution to the total mass 
of a galaxy,^ interstellar dust has a dramatic effect on the physical conditions 
and processes taking place within the universe, in particular, the evolution of 
galaxies and the formation of stars and planetary systems (see the introduction 
section of Li & Greenberg 2003). 

In this review I will concentrate on the radiative properties of interstellar 
dust. I will distinguish the dust in terms of 3 components: cold dust with 
steady-state temperatures of 15K^r^25K in thermal equilibrium with the 
solar neighbourhood interstellar radiation field, very cold dust with T<10K 
(typically T~ 5 K), and warm dust undergoing “temperature spikes” via single¬ 
photon heating. I will first summarize in §2 the observational constraints on 
the physical and chemical properties of interstellar dust. I will then discuss 
in §3 the heating and cooling properties of the warm and cold interstellar dust 
components. In §4 I will demonstrate the robustness of the silicate-graphite- 
PAHs interstellar grain model by showing that this model closely reproduces 
the observed infrared (IR) emission of the Milky Way, the Small Magellanic 
Cloud, and the ringed Sb galaxy NGC 7331. The possible existence of a pop¬ 
ulation of very cold dust (with T< 10 K) in interstellar space will be discussed 
in §5. Whether ultrasmall grains can really cool down to T <2.7 K and appear 
in absorption against the cosmic microwave background (CMB) radiation will 
be discussed in §6. 

2. Observational Constraints on Interstellar Dust 

Our knowledge of interstellar dust regarding its composition, size and shape 
is mainly derived from its interaction with electromagnetic radiation: atten¬ 
uation (absorption and scattering) and polarization of starlight, and emission 
of IR and far-IR radiation. The principal observational keys, both direct and 
indirect, used to constrain the properties of dust are the following: 
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■ (1) Interstellar Extinction . Extinction is a combined effect of absorp¬ 
tion and scattering: a grain in the line of sight between a distant star 
and the observer reduces the starlight by a combination of scattering and 
absorption (the absorbed energy is then re-radiated in the IR and far-IR). 

- The wavelength dependence of interstellar extinction - “interstel¬ 
lar extinction curve”, most commonly determined from the “pair- 
method”,^ rise from the near-IR to the near-UV, with a broad ab¬ 
sorption feature at about A“^^4.6 //m“^ (Af«2175 A), followed by 
a steep rise into the far-UV A“^f«10/im“^. —> The extinction 
curve tells us the size (and to a less extent, the composition) of 
interstellar dust. Since it is generally true that a grain absorbs and 
scatters light most effectively at wavelengths comparable to its size 
Af«27ra (Krugel 2003), there must exist in the ISM a population 
of large grains with a^A/27rRi0.1 //m to account for the extinc¬ 
tion at visible wavelengths, and a population of ultrasmall grains 
with a^A/27rKi0.016//m to account for the far-UV extinction at 
A=0.1 (see Li 2004a for details). 

- The optical/UV extinction curves show considerable regional vari¬ 
ations. —> Dust grains on different sightlines have different 
size distributions (and/or different compositions). 

- The optical/UV extinction curve in the wavelength range of 0.125 
<A<3.5 /xm can be approximated by an analytical formula involv¬ 
ing only one free parameter: Ry^Ay/E{B — V), the total-to- 
selective extinction ratio (Cardelli, Clayton, & Mathis 1989), wher¬ 
eas the near-IR extinction curve (0.9 /xm<A<3.5 /xm) can be fitted 
reasonably well by a power law A(A)~A“^'^, showing little envi¬ 
ronmental variations. 

The Galactic mean extinction curve is characterized by Ry ^ 3.1. 
Values of Ry as small as 2.1 (the high latitude translucent molec¬ 
ular cloud HD 210121; Larson, Whittet, & Hough 1996; Li & 
Greenberg 1998) and as large as 5.6 (the HD 36982 molecular 
cloud in the Orion nebula) have been observed in the Galactic re¬ 
gions. More extreme extinction curves have been reported for ex- 
tragalactic objects.^ —> The optical/UV extinction curve and 
the value of Ry depend on the environment: lower-density re¬ 
gions have a smaller Ry, a stronger 2175 A hump and a steeper 
far-UV rise (A“^>4^m“^), implying smaller grains in these re¬ 
gions; denser regions have a larger Ry, a weaker 2175 A hump 
and a flatter far-UV rise, implying larger grains. 

- In the Small Magellanic Cloud (SMC), the extinction curves of 
most sightlines display a nearly linear steep rise with A“^ and 



an extremely weak or absent 2175 A hump (Lequeux et al. 1982; 
Prevot et al. 1984). —> Grains in the SMC are smaller than 
those in the Milky Way diffuse ISM as a result of either more 
efficient dust destruction in the SMC due to its harsh environ¬ 
ment of the copious star formation associated with the SMC 
Bar or lack of growth due to the low-metallicity of the SMC, or 
both. Regional variations also exist in the SMC extinetion curves.® 

- The Large Magellanic Cloud (LMC) extinction curve is character¬ 
ized by a weaker 2175 A hump and a stronger far-UV rise than the 
average Galactic extinction curve (Nandy et al. 1981; Koornneef 
&Code 1981).^ 

(2) Interstellar Polarization . For a non-spherical grain, the light of 
distant stars is polarized as a result of differential extinction for different 
alignments of the electric vector of the radiation. 

- The interstellar polarization curve rises from the IR, has a maxi¬ 
mum somewhere in the optical and then decreases toward the UV. 
—> This tells us that (1) some fraction of interstellar grains are 
both non-spherical and aligned by some process; (2) the bulk 
of the aligned grains responsible for the peak polarization (at 
Af«0.55 /xm) have typical sizes of aRiA/27rf«0.1 ^um; and (3) the 
ultrasmall grain component responsible for the far-UV extinc¬ 
tion rise is either spherical or not aligned. 

- The optical/UV polarization curve P{X) can be closely approxi¬ 
mated by the “Serkowski law”, an empirical function: P{X)/Pma.x 
= exp[—iT ln^(A/Amax)], where the only one free parameter. A m ax , 
is the wavelength where the maximum polarization Pmax occurs, 
and K is the width parameter: FfRil.GhAmax + 0.01; A m ax is in¬ 
dicative of grain size and correlated with Ry- Ry ~ (5.6=fc0.3)A max 
(Amax is in micron; see Whittet 2003). —> The sightlines with 
larger Amax are rich in larger grains and have larger Ry for 
their extinction curves. A close correlation between Amax and 
the size of the aligned grains [e.g., Amax~ 27 ra(n — 1) for dielec¬ 
tric cylinders of radius a and refractive index n] is predicted by 
interstellar extinction calculations. 

- The near-IR (1.64;um<A<5 ^um) polarization is better approxi¬ 
mated by a power law P{X) oc X~^, with /3~1.8±0.2, independent 
of Amax (Martin & Whittet 1990, Martin et al. 1992). 

(3) Interstellar Scattering. Scattering of starlight by interstellar dust 
is revealed by reflection nebulae (dense clouds illuminated by embed¬ 
ded or nearby bright stars), dark clouds (illuminated by the general in¬ 
terstellar radiation field [ISRF]), and the diffuse Galactic light (“DGL”; 
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starlight scattered off the general diffuse ISM of the Milky Way Galaxy 
illuminated by the general ISRF). The scattering properties of dust grains 
(albedo = ratio of scattering to extinction, and phase function) provide a 
means of constraining the optical properties of the grains and are there¬ 
fore indicators of their size and composition and provide diagnostic tests 
for dust models. 

- The albedo in the near-IR and optical is quite high (~ 0.6), with a 
clear dip to ~ 0.4 around the 2175 A hump, a rise to ~ 0.8 around 
A~^f«6.6;um“^, and a drop to ~0.3 by A“^^10;um“^; the scat¬ 
tering asymmetry factor almost monotonically rises from ~ 0.6 to 
~0.8 from A“^Ril/im“^ to A“^f«10/im“^ (see Gordon 2004). 
—> An appreciable fraction of the extinction in the near-IR 
and optical must arise from scattering; the 2175 A hump is an 
absorption feature with no scattered component (see §2.4); and 
ultrasmall grains are predominantly absorptive. 

- Surprisingly high near-IR albedo has been reported for several re¬ 
gions: ~0.86 at the K' band (A?5i2.1 //m) for the prominent dust 
lane in the evil eye galaxy NGC4826 (Witt et al. 1994), ~0.9 at 
the K' band for the dust in the M51 arm (Block 1996), and ~0.7 
at the J (Af«1.26 /xm) and H (Af«1.66 /xm) bands and ~ 0.6 at the K 
(Af«2.16/xm) band for the Thumbprint Nebula (Lehtinen & Mat- 
tila 1996), in comparison with ~ 0.2 at A=2.2 ^m predicted by the 
conventional dust models (Draine & Lee 1984; Li & Greenberg 
1997). —> This implies that for NGC4826 and M51 (1) a pop¬ 
ulation of grains at least ~ 0.5 jivu in radii, twice as large as as¬ 
sumed by standard models, may exist in these environments and 
are responsible for the near-IR scattering; and/or (2) the measured 
high K' surface-brightness and the deduced high albedo may in part 
be caused by the thermal continuum emission from stochastically 
heated ultrasmall grains (Witt et al. 1994; Block 1996). For the 
Thumbprint Nebula, the high near-IR albedo is readily explained 
by grain growth (larger-than-average grain sizes) and the accretion 
of an ice mantle (see Pendleton, Tielens, & Werner 1990; §8 in Li 
& Greenberg 1997). 

- Scatterings of X-rays by interstellar dust have also been observed 
as evidenced by “X-ray halos” formed around an X-ray point source 
by small-angle scattering. The intensity and radial profile of the 
halo depends on the composition, size and morphology and the 
spatial distribution of the scattering dust particles (see Dwek et 
al. 2004 for a review). The total and differential cross sections 
for X-ray scattering approximately vary as and a®, respectively, 
the shape and intensity of X-ray halos surrounding X-ray point 



sources therefore provide one of the most sensitive constraints 
on the largest grains along the sightline, while these grains are 
gray at optieal wavelengths and therefore the near-IR to far-UV 
extinetion modeling is unahle to eonstrain their existenee. 

A reeent study of the X-ray halo around Nova Cygni 1992 by Witt, 
Smith, & Dwek (2001) pointed to the requirement of large inter¬ 
stellar grains (a~ 0.25-2 ^m), eonsistent with the reeent Ulysses 
and Galileo deteetions of interstellar dust entering our solar system 
(Grun et al. 1994; Friseh et al. 1999; Landgraf et al. 2000). But 
Draine & Tan (2003) found that the silieate-graphite-PAH model 
with the dust size distributions derived from the near-IR to far-UV 
extinetion modeling (Weingartner & Draine 2001a) and IR emis¬ 
sion modeling (Li & Draine 2001b) is able to explain the observed 
X-ray halo. 

(4) Spectroscopic Extinction and Polarization Features: The 2175 A 
Extinction Hump - the strongest spectroscopic extinction feature. 

- Its strength and width vary with environment while its peak posi¬ 
tion is quite invariant: the eentral wavelength of this feature varies 
by only ±0.46% {2a) around 2175 A (4.6 ;um“^), while its FWHM 
varies by ±12% {2a) around 469 A (f»l /rm“^). 

- Its earrier remains unidentified 39 years after its first deteetion 
(Steeher 1965). It is generally believed to be eaused by aromatie 
earbonaeeous (graphitie) materials, very likely a eosmie mixture of 
polyeyelie aromatie hydroearbon (PAH) moleeules (Joblin, Leger 
& Martin 1992; Li & Draine 2001b; Draine 2003a). 

- For most sightlines, this feature is unpolarized. So far only two 
lines of sight toward HD 147933 and HD 197770 have a weak 2175 A 
polarization feature deteeted (Clayton et al. 1992; Anderson et al. 
1996; Wolff et al. 1997; Martin, Clayton, & Wolff 1999). Even for 
these sightlines, the degree of alignment and/or polarizing ability 
of the earrier should be very small (see Li & Greenberg 2003). 

- Exeept for the deteetion of seattering in the 2175 A hump in two re- 
fleetion nebulae (Witt, Bohlin, & Steeher 1986), the 2175 A hump 
is thought to be predominantly due to absorption, suggesting its 
earrier is suffieiently small to be in the Rayleigh limit. 

- The deteetions of this feature in distant objeets have been reported 
by Malhotra (1997) in the eomposite absorption speetrum of 96 
intervening Mg II absorption systems at 0.2<2:<2.2; by Cohen et 
al. (1999) in a damped Eya absorber (DEA) at 2 = 0.94; by Toft, 
Hjorth & Burud (2000) in a lensing galaxy at 2 = 0.44; by Motta 
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et al. (2002) in a leasing galaxy at 2 ; = 0.83; and very reeently by 
Wang et al. (2004) in 3 intervening quasar absorption systems at 
1.4^z^l.5. 

■ (5) Spectroscopic Extinction and Polarization Features: The 9.7 ^um 
and 18 jivti (Silicate) Absorption Features - the strongest IR Absorp¬ 
tion features. 

- Ubiquitously seen in a wide range of astrophysieal environments, 
these features are almost eertainly due to silicate minerals: they 
are respectively ascribed to the Si-0 stretching and O-Si-0 bending 
modes in some form of silicate material (e.g. olivine Mg2a:Fe2-2a:Si04). 

- The observed interstellar silicate bands are broad and relatively 
featureless. —> Interstellar silicates are largely amorphous rat¬ 
her than crystalline. Li & Draine (2001a) estimated that the amount 
of a<l /xm crystalline silicate grains in the diffuse ISM is <5% of 
the solar Si abundance.® 

- The strength of the 9.7 jivu feature is approximately Ar 9 , 7 ^m/^V'~ 
1/18.5 in the local diffuse ISM. —> Almost all Si atoms have 
been locked up in silicate dust, if assuming solar abundance 
for the ISM.9 

- The 9.7 and IS/xm silicate absorption features are polarized in 
some interstellar regions, most of which are featureless. The only 
exception is AFGL 2591, a molecular cloud surrounding a young 
stellar object, which displays a narrow feature at 11.2 /xm superim¬ 
posed on the broad 9.7 /xm polarization band, generally attributed 
to annealed silicates (Aitken et al. 1988). 

■ (6) Spectroscopic Extinction and Polarization Features: The 3.4 fim. 
(Aliphatic Hydrocarbon) Absorption Feature. 

- Widely seen in the diffuse ISM of the Milky Way Galaxy and other 
galaxies (e.g. Seyfert galaxies and ultraluminous infrared galax¬ 
ies, see Pendleton 2004 for a recent review), this strong absorption 
band is attributed to the C-H stretching mode in aliphatic hydro¬ 
carbon dust. Its exact nature remains uncertain, despite 23 years’ 
extensive investigation with over 20 different candidates proposed 
(see Pendleton & Allamandola 2002 for a summary). So far, the 
experimental spectra of hydrogenated amorphous carbon (HAG; 
Schnaiter, Henning & Mutschke 1999, Mennella et al. 1999) and 
the organic refractory residue, synthesized from UV photoprocess¬ 
ing of interstellar ice mixtures (Greenberg et al. 1995), provide the 
best fit to both the overall feature and the positions and relative 




strengths of the 3.42 ^um, 3.48 /xm, and 3.51 subfeatures eorre- 
sponding to symmetrie and asymmetrie stretehes of C-H bonds in 
CH 2 and CH 3 groups. Pendleton & Allamandola (2002) attributed 
this feature to hydroearbons with a mixed aromatie and aliphatie 
eharaeter. 

- The 3.4 band strength for interstellar aliphatie hydroearbon 
dust is unknown. If we adopt a mass absorption eoeffieient of 
Kabs(3.4|Um)~ 1500cm^g“^ (Li & Greenberg 2002), we would 
require ~ 68 ppm C to be loeked up in this dust eomponent to 
aeeount for the loeal ISM 3.4^m feature (Ar3.4^m/^t/~l/250; 
Pendleton et al. 1994). 

- This feature is ubiquitously seen in the diffuse ISM while never 
detected in molecular clouds. Mennella et al. (2001) and Munoz 
Caro et al. (2001) argue that this ean be explained by the eom- 
petition between dehydrogenation (destruetion of C-H bonds by 
UV photons) and rehydrogenation (formation of C-H bonds by 
H atoms interaeting with the earbon dust): in diffuse elouds, re¬ 
hydrogenation prevails over dehydrogenation; in dense moleeular 
elouds, dehydrogenation prevails over rehydrogenation as a result 
of the redueed amount of H atoms and the presenee of iee mantles 
whieh inhibits the hydrogenation of the underlying earbon dust by 
H atoms while dehydrogenation ean still proeeed sinee the UV ra¬ 
diation ean penetrate the iee layers. 

- Whether the origin of the interstellar aliphatie hydroearbon dust 
oeeurs in the outflow of earbon stars or in the ISM itself is a sub- 
jeet of debate. The former gains strength from the elose similarity 
between the 3.4 |Um interstellar feature and that of a earbon-rieh 
protoplanetary nebula CRL 618 (Lequeux & Jourdain de Muizon 
1990; Chiar et al. 1998). However, the survival of the stellar-origin 
dust in the diffuse ISM is questionable (see Draine 1990). 

- So far, no polarization has been detected for this feature (Adam¬ 
son et al. 1999).^° Speetropolarimetrie measurements for both the 
9.7 ^m silieate and the 3.4 /xm hydroearbon features for the same 
sightline would allow a direet test of the silieate eore-hydroearbon 
mantle interstellar dust model (Li & Greenberg 1997), sinee this 
model prediets that the 3.4 ^m feature would be polarized if the 
9.7 |Um feature (for the same sightline) is polarized (Li & Green¬ 
berg 2002 ). 

(7) Spectroscopic Extinction and Polarization Features: The Ice Ab¬ 
sorption Features. 
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- Grains in dark molecular clouds (usually with Ay >3 mag) obtain 
ice mantles consisting of H 2 O, NH 3 , CO, CH 3 OH, CO 2 , CH 4 , 
H 2 CO and other molecules (with H 2 O as the dominant species), as 
revealed by the detection of various ice absorption features (e.g., 
H 2 O: 3.1, 6.0/im; CO: 4.67 ^m; CO 2 : 4.27, 15.2 ^m; CH 3 OH: 
3.54, 9.75/xm; NH 3 : 2.97/xm; CH 4 : 7.68/xm; H 2 CO: 5.81/xm; 
OCN”: 4.62 /xm; see Boogert & Ehrenfreund 2004 for a review). 

- Polarization has been detected in the 3.1 /xm H 2 O, the 4.67 /xm CO 
and 4.62 /xm OCN~ absorption features (e.g. see Chrysostomou et 
al. 1996). 

■ (8) The Extended Red Emission: Dust Photoluminescence. The “Ex¬ 
tended Red Emission” (ERE), widely seen in reflection nebulae, plane¬ 
tary nebulae, HII regions, the Milky Way diffuse ISM, and other galax¬ 
ies, is a far-red continuum emission in excess of what is expected from 
simple scattering of starlight by interstellar dust (see Witt & Vijh 2004). 
The ERE is characterized by a broad, featureless band between ~ 5400 A 
and 9500 A, with a width 600 A ^ FWHM ^ 1000 A and a peak of 
maximum emission at 6100 A ^ Ap ^ 8200 A, depending on the physi¬ 
cal conditions of the environment where the ERE is produced. 

- The ERE is generally attributed to photoluminescence (PE) by some 
component of interstellar dust, powered by UV/visible photons 
with a photon conversion efficiency ? 7 pl::^ 10 % (Gordon, Witt, & 
Eriedmann 1998). — > The ERE carriers are very likely in the 
nanometer size range because nanoparticles are expected to lu¬ 
minesce efficiently through the recombination of the electron-hole 
pair created upon absorption of an energetic photon, since in such 
small systems the excited electron is spatially confined and the ra¬ 
diationless transitions that are facilitated by Auger and defect re¬ 
lated recombination are reduced (see Ei 2004a). 

- The ERE carrier remains unidentified. Various candidate materi¬ 
als have been proposed, but most of them appear unable to match 
the observed ERE spectra and satisfy the high-r/PL requirement 
(Draine 2003a; Ei & Draine 2002a; Ei 2004a; Witt & Vijh 2004). 
Promising candidates include PAHs (d’Hendecourt et al. 1986) and 
silicon nanoparticles (Eedoux et al. 1998, Witt et al. 1998, Smith 
& Witt 2002), but both have their own problems (see Ei & Draine 
2002 a). 

■ (9) Spectroscopic Emission Eeatures: The 3.3, 6.2, 7.7, 8.6,11.3 /xm 
“Unidentified Infrared (UIR) Emission features. The distinctive set of 
“UIR” emission features at 3.3, 6.2, 7.7, 8 . 6 , and 11.3 /xm are seen in 





a wide variety of objeets, ineluding planetary nebulae, protoplanetary 
nebulae, refleetion nebulae, HII regions, photodissoeiation fronts, eir- 
eumstellar envelopes, and external galaxies. 

- These “UIR” emission features are now generally identified as 
vibrational modes of PAHs (Leger & Puget 1984; Allamandola, 
Tielens, & Barker 1985): C-H stretehing mode (3.3 /um), C-C 
stretehing modes (6.2 and 7.7 ;um), C-H in-plane bending mode 
(8.6/um), and C-H out-of-plane bending mode (11.3/um).^^ The 
relative strengths and precise wavelengths of these features are 
dependent on the PAH size and its ionization state whieh is eon- 
trolled by the starlight intensity, eleetron density, and gas temper¬ 
ature of the environment (Bakes & Tielens 1994, Weingartner & 
Draine 2001b, Draine & Li 2001). 

- Stochastically heated by the absorption of a single UV/visible 
photon (Draine & Li 2001; Li 2004a), PAHs, containing ~45 ppm 
C, account for ~ 20% of the total power emitted by interstellar 
dust in the Milky Way diffuse ISM (Li & Draine 2001b). 

- The excitation of PAHs does not require UV photons; long wave¬ 
length (red and far-red) photons are also able to heat PAHs to high 
temperatures so that they emit effieiently at the UIR bands. This 
is beeause the PAH electronic absorption edge shifts to longer 
wavelengths upon ionization and/or as the PAH size increases. 

- No polarization has been deteeted for the PAH emission features 
(Sellgren, Rouan, & Leger 1988). 

- The PAH absorption features at 3.3 /rm and 6.2 ^m have been 
deteeted in both loeal sourees and Galaetie Center sourees (Sehutte 
et al. 1998; Chiar et al. 2000). The strengths of these features are in 
good agreement with those predieted from the astronomieal PAH 
model (Li & Draine 2001b). 

- PAHs can be rotationally excited by a number of physieal pro- 
eesses, ineluding eollisions with neutral atoms and ions, “plasma 
drag”, and absorption and emission of photons. It is shown that 
these proeesses ean drive PAHs to rapidly rotate, with rotation 
rates reaehing tens of GHz. The rotational eleetrie dipole emission 
from these spinning PAH moleeules is eapable of aeeounting for 
the observed “anomalous” microwave emission (Draine & Lazar- 
ian 1998a,b; Draine 1999; Draine & Li 2004). 

(10) IR Emission from Interstellar Dust. Interstellar grains absorb sta¬ 
rlight in the UV/visible and re-radiate in the IR. The IR emission spec¬ 
trum of the Milky Way diffuse ISM, estimated using the IRAS 12, 25, 60 
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and 100 /xm broadband photometry, the DIRBE-COBE 2.2, 3.5, 4.9, 12, 
25, 60, 100, 140 and 240/xm broadband photometry, and the EIRAS- 
COBE 110;um<A<3000^m speetrophotometry, is eharaeterized by a 
modified blaek-body of A“^'^i?A(T=19.5 K) peaking at ~ 130 ^m in the 
wavelength range of 80^m^A^1000/um, and a substantial amount of 
emission at A^60 ^m whieh far exeeeds what would be expeeted from 
dust at Tk,20¥^ (see Draine 2003a). In addition, speetrometers aboard 
the IRTS (Onaka et al. 1996; Tanaka et al. 1996) and ISO (Manila et 
al. 1996) have shown that the diffuse ISM radiates strongly in emission 
features at 3.3, 6.2, 7.7, 8.6, and 11.3 /um. 

- The emission at A^60 ^m aeeounts for ~ 65% of the total emitted 
power. —> There must exist a population of “cold dust” in the size 
range of a>250 A, heated by starlight to equilibrium temperatures 
15K^r^25K and eooled by far-IR emission (see Ei & Draine 
2001 b). 

- The emission at A^60 nvn aeeounts for ~ 35% of the total emitted 
power. —> There must exist a population of “warm dust” in the 
size range of a <250 A, stoehastieally heated by single starlight 
photons to temperatures T»20 K and eooled by near- and mid-IR 
emission (see Ei & Draine 2001b; Ei 2004a). 

■ (11) Interstellar Depletions. Atoms looked up in dust grains are “de¬ 
pleted” from the gas phase. The dust depletion ean be determined from 
eomparing the gas-phase abundanees measured from optieal and UV ab¬ 
sorption speetroseopie lines with the assumed referenee abundanees of 
the ISM (total abundanees of atoms both in gas and in dust; also known 
as “standard abundanees”, “interstellar abundanees”, and “eosmie abun¬ 
danees”). The total interstellar abundanees are usually taken to be solar, 
although Snow & Witt (1996) argued that interstellar abundanees are 
appreeiably subsolar (~ 70% solar). Interstellar depletions allow us to 
extraet important information about the eomposition and quantity of in¬ 
terstellar dust: 

- In low density elouds. Si, Ee, Mg, C, and O are depleted. —> 
Interstellar dust must contain an appreciable amount of Si, Fe, 

Mg, C and O. Indeed, all eontemporary interstellar dust models 
eonsist of both silieates and earbonaeeous dust. 

- Erom the depletion of the major elements Si, Ee, Mg, C, and O one 
ean estimate the gas-to-dust mass ratio to be 165.13 

- In addition to the silieate dust eomponent, there must exist another 
dust population, sinee silicates alone are not able to account for 
the observed amount of extinction relative to H although Si, Mg, 
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and Fe are highly depleted in the ISM. Even if all Si, Fe, and Mg 
elements are loeked up in submieron-sized silieate grains, they ean 
only aeeount for ~ 60% of the total observed optieal extinetion.^^ 

3. Warm and Cold Dust in the ISM 

As summarized in §2, the interstellar extinetion, polarization, seattering, the 
near, mid, and far-IR emission and the 3.3-11.3/um PAH emission features 
point to the existenee of two dust populations in interstellar spaee: 

■ There exists a population of large grains with a;^250 A. Illuminated 
by the general interstellar radiation field (ISRF), these grains, defined 
as “ cold dust ”, obtain equilibrium temperatures of 15K^r^25K and 
emit strongly at wavelengths A;^60 /xm. These grains are responsible for 
the near-IR/optieal extinetion, seattering, polarization and the emission 
at A;^60 ^m. 

The equilibrium temperature T for a large grain of spherieal radius a is 
determined by balaneing absorption and emission: 

POO POO 

/ Cahs{a, X)cuxdX = / C'abs(a, A)47rHA(r)(iA , (1) 

where Cabs (a, A) is the absorption eross seetion for a grain with size a 
at wavelength A, c is the speed of light, B\{T) is the Planek funetion 
at temperature T, and u\ is the energy density of the radiation field. In 
Figure 1 we display fhese “equilibrium” femperafures for graphifie and 
silieafe grains as a funelion of size in environmenfs wifh various UV 
intensifies. 

■ There exisfs a population of ulfrasmall grains wifh a^250A. These 
grains have energy eonfenfs smaller or eomparable fo fhe energy of a 
single sfarlighf photon. As a resulf, a single phofon ean heal a very small 
grain fo a peak lemperafure mueh higher lhan ils “sleady-slale” fem- 
peralure and fhe grain Ihen rapidly eools down fo a lemperafure below 
fhe “sleady-slale” lemperafure before Ihe nexl photon absorption event 
Sloehaslie heating by absorption of slarlighl Iherefore resulls in Iransienl 
“lemperafure spikes”, during whieh mueh of Ihe energy deposited by Ihe 
slarlighl photon is reradialed in Ihe IR - beeause of Ibis, we eall Ihese 
ulfrasmall grains “ warm dust ”.^^ These grains are responsible for Ihe 
far-UV extinetion rise and Ihe emission al A^60 iivu (ineluding Ihe 3.3- 
11.3 /xm PAH emission fealures), dominate Ihe pholoeleelrie heating of 
interstellar gas (see §2.5 in Fi 2004a), and provide mosl of Ihe grain 
surfaee area in Ihe diffuse ISM. 

Sinee ulfrasmall grains will nol mainlain “equilibrium lemperalures”, 
we need to ealeulale Iheir lemperalure (energy) probabilily dislribulion 
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Figure 1. Equilibrium temperatures for graphite (dotted lines) and silicate grains (solid lines) in envi¬ 
ronments with various starlight intensities (in units of the Mathis, Mezger, & Panagia 1983 [MMP] solar 
neighbourhood ISRF). Taken from Li & Draine (2001b). 


functions in order to calculate their time-averaged IR emission spectrum. 
There have been a number of studies on this topic since the pioneering 
work of Greenberg (1968). A recent extensive investigation was carried 
out by Draine & Li (2001). We will not go into details in this review, but 
just refer those who are interested to Draine & Li (2001) and a recent 
review article of Li (2004a). 

For illustration, we show in Figure 2 the energy probability distribu¬ 
tion functions dP/d\nE (where dP is the probability that a grain will 
have vibrational energy in interval [E,E + dE]) for PAHs with radii 
a = 5,10, 25, 50, 75,100,150, 200, 300A illuminated by the general 
ISRF. It is seen that very small grains (a^lOOA) have a very broad 
P{E), and the smallest grains (a<30 A) have an appreciable probability 
Pq of being found in the vibrational ground state E=0. As the grain size 
increases, P{E) becomes narrower, so that it can be approximated by a 
delta function for a>250 A.^® However, for radii as large as a=200A, 
grains have energy distribution functions which are broad enough that 
the emission spectrum deviates noticeably from the emission spectrum 
for grains at a single “steady-state” temperature T, as shown in Fig- 
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Figure 2. The energy probability distribution functions for charged carbonaceous grains (a = 5 A 
[C 60 H 20 ], 10 A [C 480 H 120 ], 25 A [C 7200 H 1800 ], 50, 75, 100, 150, 200, 250, 300 A) illuminated by the 
general ISRF. The discontinuity in the 5, 10, and 25 A curves is due to the change of the estimate for grain 
vibrational “temperature” at the 20th vibrational mode (see Draine & Li 2001). For 5, 10, and 25A a dot 
indicates the first excited state, and Pq is the probability of being in the ground state. Taken from Li & 
Draine (2001b). 


ure 3. For accurate computation of IR emission spectra it is therefore 
important to properly calculate the energy distribution function P{E), 
including grain sizes which are large enough that the average thermal 
energy content exceeds a few eV. 

■ Nano-sized interstellar diamond and TiC grains have been identified in 
primitive meteorites based on isotopic anomaly analysis. Exposed to 
the general ISRF, these grains are subject to single-photon heating and 
should of course be classified as “warm dusf”. Buf we should nofe fhaf 
fhey are nol represenfafive of fhe bulk infersfellar dusf (see §5.4 and §5.5 
in Li 2004a). The carriers of fhe ERE (see §2.8) and fhe 2175 A ex- 
fincfion hump (see §2.4) are also in fhe single-phofon heating regime. 
Therefore, fhey can also be classified as “warm dusf”. As a maffer of 
facf, fhe latter is aflribuled fo PAHs (see §2.4). On fhe ofher hand, fhe 
Ulysses and Galileo spacecraffs have defecfed a subsfanfial number of 
large infersfellar grains wifh a^l ^m (Grun ef al. 1994), much higher 
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Figure 3. Infrared emission spectra for small carbonaceous grains of various sizes heated by the gen¬ 
eral ISRF, calculated using the full energy distribution function P{E) (solid lines); also shown (broken 
lines) are spectra computed for grains at the “equilibrium” temperature T. Transient heating effects lead to 
significantly more short wavelength emission for a<200 A. Taken from Li & Draine (2001b). 


than expected for the average interstellar dust distribution. In the ISM, 
these should of course be considered as “cold dust” or “very cold dust” 
with r<10K. But we note that the reported mass of these large grains 
is difficult to reconcile with the interstellar extinction and interstellar el¬ 
emental abundances. 

4. The Silicate-Graphite-PAHs Interstellar Dust Model 

Various models have been proposed for interstellar dust (see Li & Greenberg 
2003, Li 2004a, Draine 2004 for recent reviews). In general, these models fall 
into three broad categories: the silicate core-carbonaceous mantle model (Li 
& Greenberg 1997), the silicate-graphite-PAHs model (Li & Draine 2001b, 
Weingartner & Draine 2001a) and the composite model (Mathis 1996, Zubko, 
Dwek & Arendt 2004). In this review I will focus on the IR emission calculated 
from the silicate-graphite-PAHs model and refer those who are interested in a 
detailed comparison between different models to my recent review articles (Li 
2004a, Li & Greenberg 2003). 
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The silicate-graphite-PAHs model, consisting of a mixture of amorphous 
silicate dust and carhonaceous dust - each with an extended size distribution 
ranging from molecules containing tens of atoms to large grains ;^1 /rm in di¬ 
ameter, is a natural extension of the classical silicate-graphite model (Mathis, 
Rumpl, & Nordsieck 1977; Draine & Lee 1984). We assume that the carhona¬ 
ceous grain population extends from grains with graphitic properties at radii 
a^50 A, down to particles with PAH-like properties at very small sizes. 

With the temperature (energy) prohahility distribution functions calculated 
for ultrasmall grains undergoing “temperature spikes” and equilibrium tem¬ 
peratures calculated for large grains illuminated by the local ISRF, the silicate- 
graphite-PAHs model with grain size distributions consistent with the observed 
Ry=3.1 interstellar extinction (Weingartner & Draine 2001a), is able to repro¬ 
duce the observed near-IR to submillimeter emission spectrum of the diffuse 
ISM, including the PAH emission features at 3.3, 6.2, 7.7, 8.6, and 11.3 ;um. 
This is demonstrated in Figure 4 and Figure 5 for the high-latitude “cirrus” 
cloud and 2 regions in the Galactic plane (see Li & Draine 2001b for details). 



Figure 4. Comparison of the model to the observed emission from the diffuse ISM at high galactic 
latitudes (|6|>25°). Curves labelled Bgii and Bcarb show emission from “big” (a>250A) silicate and 
carbonaceous grains; curves labelled Sgii and Scarb show emission from “small” (a<250 A) silicate and 
carbonaceous grains (including PAHs). Triangles show the model spectrum (solid curve) convolved with 
the DIRBE filters. Observational data are from DIRBE (diamonds) and FIRAS (squares). Taken from Li & 
Draine (2001b). 
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Figure 5. Infrared emission from dust plus starlight for two regions in the Galactic plane: (a) the MIRS 
region (44° < I < 44°40', -0°40' <b< 0°), and (b) the MRS region (47°30' < I < 48°, \b\ < 15')- 
The starlight intensity heating the dust has been taken to be twice the MMP ISRF. The solid curve shows 
the overall model spectrum; triangles show the model spectrum convolved with the DIRBE filters. DIRBE 
observations are shown as diamonds. For the MIRS field we show the IRTS MIRS 5-12 fim spectrum (thin 
solid line). For the MRS field we show the IRTS NIRS 2.8-3.9 fim spectrum (thin solid line, also shown as 
cross-dotted curve in inset). Taken from Li & Draine (2001b). 


The silicate-graphite-PAHs model, with size distributions consistent with 
the SMC Bar extinction curve (Weingartner & Draine 2001a), is also successful 
in reproducing the observed IR emission from the SMC (Li & Draine 2002c), 
as shown in Figure 6. The dust in the SMC is taken to be illuminated by a 
distribution of starlight intensities. Following Dale et al. (2001), we adopt a 
simple power-law function for the starlight intensity distribution. The SMC, 
with a low metallicity (~ 10% of solar) and a low dust-to-gas ratio (~ 10% of 
the Milky Way), has a very weak or no 2175 A extinction hump in its extinc¬ 
tion curves for most sightlines (see §2.1) and very weak 12 fim emission (see 
Fig. 6) which is generally attributed to PAHs, supporting the idea of PAHs as 
the carrier for the 2175 A extinction hump.^^ 
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Figure 6. Comparison of the model (solid line) to the observed emission from the SMC obtained 
by COBE/DIRBE (diamonds) and IRAS (squares) averaged over a 6.25 deg^ region including the optical 
bar and the Eastern Wing. Triangles show the model spectrum convolved with the DIRBE filters. Stellar 
radiation (dot-dashed line) dominates for A<6 fim. Grains are illuminated by a range of radiation intensities 
dNn/dU oc with7V^°* 5.4 x 10^^ cm-^. Taken from Li & Draine (2002c). 


Very recently, the silicate-graphite-PAHs model has also been successfully 
applied to NGC7331, a ringed Sb galaxy. Using the same set of dust pa¬ 
rameters determined for the Milky Way diffuse ISM (i?v/=3.1; Weingartner & 
Draine 2001a), as shown in Figure 7, this model fits the IR emission observed 
by the IRAC instrument at 3.6, 4.5, 5.8 and 8/xm and the MIPS instrument 
at 24, 70 and 160^m aboard the Spitzer Space Telescope and the 450 and 
850 /xm SCUBA submillimeter emission observed by JCMT, both for the ring 
and inside star-forming region and for the galaxy as whole (see Regan et al. 
2004 for details). The model also closely reproduces the observed 6.2, 7.7, 
8 .6, 11.3 and 12.7 fj,m. PAH emission features (see Fig. 2 of Smith et al. 2004). 

5. Very Cold Dust or How Cold Could Galaxies Be? 

In the 1996 South African “Cold Dust” Symposium (Block & Greenberg 
1996), the possible existence of a population of very cold dust (with equilib¬ 
rium temperatures T<10K) in interstellar space was a subject received much 
attention. In his invited paper titled “How Cold Could Galaxies Be?” published 
in the proceedings for that symposium, Mike Disney wrote “... An eminent cos- 
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Figure 7. IR emission and model fits to the NGC7331 ring (a) and the entire galaxy (b). The thick 
solid lines and triangles are the model-predicted fluxes, and the squares are the observed fluxes. The broken 
curves indicate the contributions of the different model components. Taken from Regan et al. (2004). 


mologist once advised me to forget all about very cold dust because the law 
ensures that it cannot emit, and therefore by implication cannot absorb, much 
radiation. He sounded plausible, as Cosmologists are apt to sound, but he was 
in fact totally wrong, as Cosmologists are apt to be.” 

Disney (1996) argued that for grains with higher IR and far-IR emissivities, 
they can achieve rather low temperatures. At first glance, this appears plausible 
as can be seen in Eq.(l): for a given interstellar radiation field, grains wifh 
fixed UV and opfical absorpfion properties would obfain lower femperafures 
if fheir long wavelengfh emissivities are enhanced. Therefore, Disney (1996) 
wrofe “... [Since] we are not confident about their [interstellar grains] size 
distribution and their emissivities, particularly at long wavelength, ... we have 
to keep our minds open to the possible existence of a significant amounts 
of very cold (<10K) material in spiral galaxies.” 
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So far, the detection of very cold dust (r<10K) in interstellar space on 
galactic scales has been reported for various objects: NGC4631, a low metal- 
licity (~ 1/2 of solar) interacting galaxy with T ~ 4-6 K (Dumke, Krause, & 
Wielebinski 2004); NGC 1569, a low metallicity (~ 1/4 of solar) starburst¬ 
ing dwarf galaxy with T~5-7K (Galliano et al. 2003); inactive spiral galax¬ 
ies UGC3490 with r~9K (Chini et al. 1995), NGC 6156 with T~8.6K, 
and NGC 6918 with T~9.4K (Krugel et al. 1998); and several irregular and 
blue compact dwarf galaxies with T<10 K in the Virgo Cluster (Popescu et al. 
2002). Very cold dust with T~ 4-7 K has also been detected in the Milky Way 
Galaxy (Reach et al. 1995; also see Boulanger et al. 2002). This component 
is widespread and spatially correlated with the warm component (16-21 K). 
By comparing the dust mass calculated from the IRAS data with the molecular 
and atomic gas masses of 58 spiral galaxies, Devereux & Young (1990) argued 
that the bulk of the dust in spiral galaxies is <15 K regardless of the phase of 
the ISM. 

How can dust get so cold? In literature, suggested solutions include (1) the 
dust is deeply embedded in clumpy clouds and heated by the far-IR emission 
from “classical grains” (Galliano et al. 2003; Dumke et al. 2004);^® (2) the 
dust has unusual optical properties (e.g. fractal or porous grains with enhanced 
submm and mm emissivity; Reach et al. 1995; Dumke et al. 2004). However, 
as discussed in detail by Li (2004b), while the former solution appears to 
be inconsistent with the fact that the very cold dust is observed on galactic 
scales, the latter violates the Kramers-Kronig dispersion relation (Purcell 
1969; Draine 2003b; Li 2003b), except for extremely elongated conducting 
dust (Li 2003a). Perhaps the submm and mm excess emission (usually at¬ 
tributed to very cold dust) is from something else? To avoid the “temperature” 
problem, one can adopt the Block direct method to identify and characterize 
the presence and distribution of the cold and very cold dust: using near-IR 
camera arrays and subtracting these from optical CCD images. This method 
measures the dust extinction cross section and does not require the knowledge 
of dust temperature (see Block 1996; Block et al. 1994a,b, 1999). Otherwise, a 
detailed radiative transfer treatment of the interaction of the dust with starlight 
(e.g. Popescu et al. 2000, Tuffs et al. 2004) together with a physical interstellar 
dust model (e.g. Li & Draine 2001b,2002c) is required. 

6. Can Dust Get Down to 2.7 K and Appear in Absorption against 
the CMB? 

Ultrasmall grains spend most of their time at their vibrational ground state dur¬ 
ing the interval of two photon absorption events (§3 and Fig. 2; Draine & Li 
2001). In the 1996 South African ‘"Cold Dust” Symposium, it was heatedly ar¬ 
gued that these grains could obtain a vibrational temperature less than the 2.7 K 
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temperature of the CMB so that they eould be deteeted in absorption against 
the CMB (Duley & Poole 1998). However, based on detailed modeling of the 
exeitation and de-exeitation of these grains, we found that even though these 
grains do have a large population in the vibrational ground state, nevertheless 
the vibrational levels are suffieiently exeited that the grains would appear in 
emission against the CMB with brightness temperature ^9 K (see Draine & Li 
2004b for details). 
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Notes 


1. Struve & Elvey (1938) appear to be first to show that the interstellar gas is mostly hydrogen. 

2. Let interstellar grains be approximated by a single size of a (spherical radius) with a column density 
of Nfi. The gas-to-dust mass ratio is 

rugas _ (PH + [He/H]ig^^ ^He) Ah ^ 1.4/rH Ah qx 

mdust (4/3) Tva^pdNd (4/3) pdNd 

where fin and fine are respectively the atomic weight of H and He; [He/H]jgj^ is the interstellar He 
abundance (relative to H), which is taken to be that of the solar value, ~0.1; pu is the mass density of 
interstellar dust. The hydrogen (of all forms)-to-dust column density Nn/N^ can be derived from 


_ Nn _ _ Ry Nn _ 

E{B-V) ~ Av/Rv ~ 1.0867ra2Qext(V)A^d 


( 3 ) 


where Ry is the total-to-selective extinction ratio, Qext(V^) is the dust extinction efficiency at V-band 
(A = 5500 A). Therefore, the gas-to-dust ratio can be readily estimated from 


22s^ ~ 1 14 Qext(L) 1 Ah .X. 212 / Qext[V] \ / 0.1/jm \ / 2.5gcm 

mdust apd Rv - V) V 1.5 / V a / ( Pd ) 

( 4 ) 

If we take canonical numbers of /?v~3.1, Qext(V^)~l-5, a=0.1 pm, pd~2.5gcm mgas/mdust 
would be around 210. 

3. In our Milky Way Galaxy, interstellar matter (gas and dust; 7 di 3 x IO^Mq), contributes roughly 
20% of the total stellar mass (4 ± 2 x IO^^Mq). Therefore, the mass fraction of interstellar dust is just 

~0.1% in our Galaxy ('-^ 1 di 0.2 x IO^^Mq within lOkpc, Kennicutt 2001)! 

4. In this method, the wavelength dependence of interstellar extinction is obtained by comparing the 
spectra of two stars of the same spectral type, one of which is reddened and the other unreddened. 

5. Falco et al. (1999) found Ry^lA for an elliptical lensing galaxy at 2:^:^0.96, and Ry^ 7.2 for a 
spiral lensing galaxy at 2 /^0.68. Wang et al. (2004) found the extinction curves for two intervening quasar 
absorption systems at z^l.5 to have Ry^0.1,\.9. 

6 . For example, there is at least one line of sight (Sk 143=AvZ 456) with an extinction curve with a 
strong 2175A hump detected (Lequeux et al. 1982; Prevot et al. 1984; Bouchet et al. 1985; Thompson et al. 
1988; Gordon & Clayton 1998). This sightline passes through the SMC wing, a region with much weaker 
star formation (Gordon & Clayton 1998). The sight lines which show no 2175A hump all pass through the 
SMC Bar regions of active star formation (Prevot et al. 1984; Gordon & Clayton 1998). 
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7. Strong regional variations in extinction properties have also been found in the LMC (Clayton & 
Martin 1985; Fitzpatrick 1985,1986; Misselt, Clayton, & Gordon 1999): the sightlines toward the stars 
inside or near the supergiant shell, LMC 2, which lies on the southeast side of the 30 Dor star-forming 
region, have very weak 2175 A hump (Misselt et al. 1999). 

8. Kemper, Vriend & Tielens (2004) found that crystalline fraction of the interstellar silicates along the 
sightline towards the Galactic Center is ~ 0.2%. 

9. The amount of Si (relative to H) required to deplete in dust to account for the observed 9.7 
feature strength is 


.H. 


Arg.T^m 1 _ Arg.r^m Av 1 

A^h Kj^®(9.7/2m)/2sii Ay Nh fim) fisil 


(5) 


where /cj^®(9.7/im) is the silicate mass absorption coefficient at A=9.7 /im; ^gil is the silicate molecular 
weight. With «g^®(9.7 /im):=::2850 cm^ and /2.sii~172/iH for amorphous olivine MgFeSi 04 , the local 
diffuse ISM (Arg.y ^ni/Av^l/18.5, Ay /Ah~5.3x 10“^^ mag cm^) requires 1^35 ppm. 

10. Hough et al. (1996) reported the detection of a weak 3.47 / 2 .m polarization feature in the Becklin- 
Neugebauer object in the OMC-1 Orion dense molecular cloud, attributed to carbonaceous materials with 
diamond-like structure. See Li (2004a) and Jones & d’Hendecourt (2004) for a detailed discussion on 
interstellar diamond. 

11. Other C-H out-of-plane bending modes at 11.9, 12.7 and 13.6/ 2 .m have also been detected. The 
wavelengths of the C-H out-of-plane bending modes depend on the number of neighboring H atoms: 

11.3 /im for solo-CH (no adjacent H atom), 11.9 /im for duet-CH (2 adjacent H atoms), 12.7 /im for trio- 
CH (3 adjacent H atoms), and 13.6 /im for quartet-CH (4 adjacent H atoms). 

12. Li & Draine (2002b) have modeled the excitation of PAH molecules in UV-poor regions. It was 
shown that the astronomical PAH model provides a satisfactory fit to the UIR spectrum of vdB 133, a 
reflection nebulae with the lowest ratio of UV to total radiation among reflection nebulae with detected UIR 
band emission (Uchida, Sellgren, & Wemer 1998). 

13. Let [X/H]q be the interstellar abundance of Xrelative to H (we assume interstellar abundances to be 

those of the solar values: [C/H]0p=:391 parts per million [ppm], [O/H]q^ 501 ppm, [Mg/H]Q^34.5 ppm, 
[Fe/H] 0^34.4 ppm, and [Si/H]Q?^28.1 ppm [Sofia 2004]); [X/H]g^g be the amount of X in gas phase 
([C/H]gg^g^l30 ppm, [0/H]gg^g?^;375 ppm; Fe, Mg and Si are highly depleted in dust: [Fe/H]gg^g^l ppm, 
[Mg/H]gg^g^2ppm, and [Si/H]g^gP:;2 ppm [Sofia 2004]); [X/Hj^^^^ be the amount of X contained in 
dust ([C/H]j^^t = [C/H]Q - [C/H]g^^si261ppm, [0/H]j^^t~126ppm, ^32.5ppm, 

[Fe/H]^^g^:=::27.1 ppm, [Si/H]^^g^^32.4 ppm). Assuming H/C=0.5 for intersteUar carbon dust, the gas- 
to-dust mass ratio is 

tttdust P/^ldust 

where the summation is over Si, Mg, Fe, C, O and H, and fix is the atomic weight of X in unit of /rjj ~ 
1.66 X 10-24 g. 

14. Assuming all Si, Mg, and Fe elements of solar abundances are condensed in silicate dust of a sto¬ 
ichiometric composition MgFeSi 04 with a characteristic size a^O.l /im, the contribution of the silicate 
dust to the optical extinction is 



1.086 TTa^Qe^tiV)Ngil/Nii 

1.086 7ra2Qext(F) (l]x=Si,Mg,Fe X/HIqMx + 4 [Si/Hlgpo) m 
(4/3) TTa^psii 

3.2 X magcm^ . (7) 


where is the column density of silicate dust, psii^3.5gcm^ is the mass density of silicate mate¬ 
rial, and QextiV) is the visual extinction efficiency of submicron-sized silicate dust which is taken to be 
gext(U)P^1.5. 
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15. The idea of transient heating of very small grains was first introduced by Greenberg (1968). This 
process was not observed until many years later when the detection of the near-IR emission of reflection 
nebulae (Sellgren, Werner, & Dinerstein 1983) and detection by IRAS of 12 and 25 //m Galactic emission 
(Boulanger & Perault 1988) were reported. 

16. This is because for large grains individual photon absorption events occur relatively frequently and 
the grain energy content is large enough that the temperature increases induced by individual photon ab¬ 
sorptions are relatively small. 

17. Li & Draine (2002c) placed an upper limit of ~0.4% of the SMC C abundance on the amount 
of PAHs in the SMC Bar. But we note that the PAH emission features have been seen in SMC Bl#l, a 
quiescent molecular cloud (Reach et al. 2000). For this region, Li & Draine (2002c) estimated that '^3% of 
the SMC C abundance to be incorporated into PAHs. 

18. Siebenmorgen et al. (1999) argued that the dust embedded in UV-attenuated clouds within the optical 
disk of typical inactive spiral galaxies cannot become colder than ^ 6 K. 

References 

Adamson, A.J., Whittet, D.C.B., Chrysostomou, A., Hough, J.H., Aitken, D.K., Wright, G.S., 
& Roche, RF. 1999, ApJ, 512, 224 

Aitken, D.K., Roche, P.F., Smith, C.H., James, S.D., & Hough, J.H. 1988, MNRAS, 230, 629 
Allamandola, L.J., & Hudgins, D.M. 2003, in Solid State Astrochemistry, ed. V. Pirronello, J. 

Krelowski, & G. Manica (Dordrecht: Kluwer), 251 
Allamandola, LJ., Tielens, A.G.G.M., & Barker, J.R. 1985, ApJ, 290, L25 
Allamandola, L.J., Hudgins, D.M., & Sandford, S.A. 1999, ApJ, 511, L115 
Andersen, A.C., Sotelo, J.A., Niklasson, G.A., & Pustovit, V.N. 2004, see WCD04, 709 
Anderson, C.M., et al. 1996, AJ, 112, 2726 
Bakes, E.L.O., & Tielens, A.G.G.M. 1994, ApJ, 427, 822 
Barnard, E.E. 1919, ApJ, 49, 1 
Bendo, G.J., et al. 2003, AJ, 125, 2361 
Block, D.L. 1996, see BG96, 1 

Block, D.L., & Greenberg, J.M., ed., 1996, New Extragalactic Perspectives in the New South 
Africa (Dordrecht: Kluwer) (hereafter BG96) 

Block, D.L., Puerari, 1., Frogel, J.A., Eskridge, P.B., Stockton, A., & Fuchs, B. 1999, Ap&SS, 
269, 5 

Block, D.L., Berlin, G., Stockton, A., Grosbol, P., Moorwood, A.F.M., & Peletier, R.F. 1994a, 
A&A, 288, 365 

Block, D.L., Witt, A.N., Grosbol, P, Stockton, A., & Moneti, A. 1994b, A&A, 288, 383 
Bohlin, R.C., Savage, B.D., & Drake, J.R 1978, ApJ, 224, 132 
Boogert, A.C.A., & Ehrenfreund, P. 2004, see WCD04, 547 

Bouchet, R, Lequeux, J., Maurice, E., Prevot, L., & Prevot-Burnichon, M. L. 1985, A&A, 149, 
330 

Boulanger, F, & Perault, M. 1988, ApJ, 330, 964 

Boulanger, F, Bourdin, H., Bernard, J.P, & Lagache, G. 2002, in EAS Publ. Sen, Vol. 4, Infrared 
and Submillimeter Space Astronomy, ed. M. Giard, J.P. Bernard, A. Klotz, & 1. Ristorcelli 
(Paris: EDP Sciences), 151 

Bowey, J.E., & Adamson, A.J. 2002, MNRAS, 334, 94 

Cardelli, J.A., Clayton, G.C., & Mathis, J.S. 1989, ApJ, 345, 245 

Chiar, J.E., Pendleton, Y.J., Geballe, T.R., & Tielens, A.G.G.M. 1998, ApJ, 507, 281 

Chiar, J.E., et al. 2000, ApJ, 537, 749 

Chini, R., Krugel, E., Lemke, R., & Ward-Thompson, D. 1995, A&A, 295, 317 
Chrysostomou, A., Hough, J.H., Whittet, D.C.B., Aitken, D.K., Roche, P. F, & Lazarian, A. 
1996, ApJ, 465, L61 



24 


Clayton, G.C., & Martin, P.G. 1985, ApJ, 288, 558 
Clayton, G.C., et al. 1992, ApJ, 385, L53 

Cohen, R.D., Burbidge, E.M., Junkkarinen, V.T., Lyons, R.W., & Madejski, G. 1999, BAAS, 
31, 942 

Dale, D.A., Helou, G., Contursi, A., Silbermann, N.A., c& Kolhatkar, S. 2001, ApJ, 549, 215 
Devereux, N.A., & Young, J.S. 1990, ApJ, 359, 42 
Disney, M. 1996, see BG96, 21 

d’Hendecourt, L.B., Leger, A., Olofson, G., & Schmidt, W. 1986, A&A, 170, 91 
Douglas, A.E., & Herzberg, G. 1941, ApJ, 94, 381 

Draine, B.T. 1990, in ASP Conf. Ser. 12, The Evolution of the Interstellar Medium, ed. L. Blitz 
(San Erancisco: ASP), 193 

Draine, B.T. 1999, in Proc. of the EC-TMR Conf. on 3K Cosmology, ed. L. Maiani, F. Mel- 
chiorri, & N. Vittorio (Woodbury: AIP), 283 
Draine, B.T. 2003a, ARA&A, 41, 241 

Draine, B.T. 2003b, in The Cold Universe, Saas-Eee Advanced Course Vol. 32, ed. D. Pfenniger 
(Berlin: Springer-Verlag), 213 
Draine, B.T. 2004, see WCD04, 691 
Draine, B.T., & Lee, H.M. 1984, ApJ, 285, 89 
Draine, B.T., & Lazarian, A. 1998a, ApJ, 494, L19 
Draine, B.T., & Lazarian, A. 1998b, ApJ, 508, 157 
Draine, B.T., & Li, A. 2001, ApJ, 551, 807 
Draine, B.T., & Li, A. 2004a, in preparation 
Draine, B.T., & Li, A. 2004b, in preparation 
Draine, B.T., & Tan, J.C. 2003, ApJ, 594, 347 
Duley, W.W., & Poole, G. 1998, ApJ, 504, LI 13 
Dumke, M., Krause, M., & Wielebinski, R. 2004, A&A, 414, 475 
Dwek, E., Zubko, V., Arendt, R.G., & Smith, R.K. 2004, see WCD04, 499 
Falco, E.E., et al. 1999, ApJ, 523, 617 
Fitzpatrick, E.L. 1985, ApJ, 299, 219 
Fitzpatrick, E.L. 1986, AJ, 92, 1068 
Frisch, PC., et al. 1999, ApJ, 525, 492 

Galliano, F, Madden, S.C., Jones, A.P, Wilson, C.D., Bernard, J.-P, Le Peintre, F. A&A, 407, 
159 

Gordon, K.D. 2004, see WCD04, 77 

Gordon, K.D., & Clayton, G.C. 1998, ApJ, 500, 816 

Gordon, K.D., Witt, A.N., & Friedmann, B.C. 1998, ApJ, 498, 522 

Greenberg, J.M. 1968, in Stars and Stellar Systems, Vol. VII, ed. B.M. Middlehurst, & L.H. 
Aller, (Chicago: Univ. of Chicago Press), 221 

Greenberg, J.M., Li, A., Mendoza-Gemez, C.X., Schutte, W.A., Gerakines, P.A., & de Groot, 
M. 1995, ApJ, 455, L177 

Grun, E., Gustafson, B.A.S., Mann, L, et al. 1994, A&A, 286, 915 
Hartmann, J. 1904, ApJ, 19, 268 

Henning, Th., Jager, C., & Mutschke, H. 2004, see WCD04, 603 

Hough, J.H., Chrysostomou, A., Messinger, D.W., Whittet, D.C.B., Aitken, D.K., & Roche, PF 
1996, ApJ, 461, 902 

Joblin, C., Leger, A., & Martin, P. 1992, ApJ, 393, L79 
Jones, A.P, & d’Hendecourt, L.B. 2004, see WCD04, 589 
Kemper, F, Vriend, W.J., & Helens, A.G.G.M. 2004, ApJ, 609, 826 

Kennicutt, R.C. 2001, in Tetons 4: Galactic Structure, Stars and the Interstellar Medium, ed. 
C.E. Woodward, M.D. Bicay, & J.M. Shull (San Francisco: ASP), 2 



The Warm, Cold and Very Cold Dusty Universe 


25 


Koornneef, J., & Code, A.D. 1981, ApJ, 247, 860 

Kragel, E. 2003, The Physics of Interstellar Dust (Bristol: loP) 

Krugel, E., Siehenmorgen, R., Zota, V., & Chini, R. 1998, A&A, 331, L9 
Landgraf, M., Baggaley, W.J., Grun, E., Kruger, H., & Linker!, G. 2000, J. Geophys. Res., 105, 
10343 

Larson, K.A., Whittet, D.C.B., &. Hough, J.H. 1996, ApJ, 472, 755 

Ledoux, G., et al. 1998, A&A, 333, L39 

Leger, A., & Puget, J.L. 1984, A&A, 137, L5 

Lehtinen, K., & Mattila, K. 1996, A&A, 309, 570 

Lequeux, J., & Jourdain de Muizon, M. 1990, A&A, 240, L19 

Lequeux, J., Maurice, E., Prevot-Bumichon, M.-L., Prevot, L., & Rocca-Volmerange, B. 1982, 
A&A, 113, L15 
Li, A. 2003a, ApJ, 584, 593 
Li, A. 2003b, ApJ, 599, L45 

Li, A. 2004a, in ASP Conf. Ser. 309, Astrophysics of Dust, ed. A.N. Witt, G.C. Clayton, & B.T. 

Draine (San Erancisco: ASP), 417 
Li, A. 2004b, to be submitted to ApJ 
Li, A., & Draine, B.T. 2001a, ApJ, 550, L213 
Li, A., & Draine, B.T. 2001b, ApJ, 554, 778 
Li, A., & Draine, B.T. 2002a, ApJ, 564, 803 
Li, A., & Draine, B.T. 2002b, ApJ, 572, 232 
Li, A., & Draine, B.T. 2002c, ApJ, 576, 762 
Li, A., & Greenberg, J.M. 1997, A&A, 323, 566 
Li, A., & Greenberg, J.M. 1998, A&A, 339, 591 
Li, A., & Greenberg, J.M. 2002, ApJ, 577, 789 

Li, A., & Greenberg, J.M. 2003, in Solid State Astrochemistry, ed. V. Pirronello, J. Krelowski, 
& G. Manice (Dordrecht: Kluwer), 37 
Malhotra, S. 1997, ApJ, 488, LOl 
Martin, P.G., & Whittet, D.C.B. 1990, ApJ, 357, 113 
Martin, P.G., Clayton, G.C., & Wolff, M.J. 1999, ApJ, 510, 905 
Martin, P.G., et al. 1992, ApJ, 392, 691 
Mathis, J.S. 1996, ApJ, 472, 643 

Mathis, J.S., Mezger, P.G., & Panagia, N. 1983, A&A, 128, 212 
Mathis, J.S., Rumpl, W., & Nordsieck, K.H. 1977, ApJ, 217, 425 
Mattila, K., et al. 1996, A&A, 315, L353 
McKellar, A. 1940, PASP, 52, 187 

Mennella, V., Brucato, J.R., Colangeli, L., & Palumbo, P. 1999, ApJ, 524, L71 
Mennella, V., et al. 2001, A&A, 367, 355 

Misselt, K.A., Clayton, G.C., & Gordon, K.D. 1998, ApJ, 515, 128 
Motta, V., et al. 2002, ApJ, 574, 719 

Munoz Caro, G.M., Ruiterkam, R., Schulte, W.A., Greenberg, J.M., & Mennella, V. 2001, A&A, 
367, 347 

Nandy, K., Morgan, D.H., Willis, A.J., Wilson, R., & Gondhalekar, P. M. 1981, MNRAS, 196, 
955 

Onaka, T., et al. 1996, PASJ, 48, L59 
Pendleton, Y.J. 2004, see WCD04, 573 
Pendleton, Y.J., & Allamandola, L.J. 2002, ApJS, 138, 75 
Pendleton, Y.J., Tielens, A.G.G.M., & Werner, M.W. 1990, ApJ, 349, 107 
Pendleton, Y.J., Sandford, S.A., Allamandola, L.J., Tielens, A.G.G.M., & Sellgren, K. 1994, 
ApJ, 437, 683 




26 


Popescu, C.C., Misiriotis, A., Kylafis, N.D., Tuffs, R.J., Fischera, J. 2000, A&A, 362, 138 
Popescu, C.C., Tuffs, RJ., Velk, HJ., Pierini, D., & Madore, B.F. 2002, ApJ, 567, 221 
Prevot, M.L., Lequeux, J., Prevot, L., Maurice, E., & Rocca-Volmerange, B. 1984, A&A, 132, 
389 

Purcell, E.M. 1969, ApJ, 158, 433 

Reach, W.T, Boulanger, E, Contursi, A., & Lequeux, J. 2000, A&A, 361, 895 

Reach, W.T, et al. 1995, ApJ, 451, 188 

Regan, M.W., et al. 2004, ApJS, 154, 204 

Schnaiter, M., Henning, Th., & Mutschke, H. 1999, ApJ, 519, 687 

Schutte, W.A., et al. 1998, A&A, 337, 261 

Sellgren, K., Rouan, D., & Leger, A. 1988, A&A, 196, 252 

Sellgren, K., Werner, M.W., & Dinerstein, H.L. 1983, ApJ, 271, L13 

Siebenmorgen, R., Kmgel, E., & Chini, R. 1999, A&A, 351, 495 

Smith, J.D., et al. 2004, ApJS, 154, 199 

Smith, T.L., & Witt, A.N. 2002, ApJ, 565, 304 

Snow, TP, & Witt, A.N. 1996, ApJ, 468, L65 

Sofia, U.J. 2004, see WCD04, 393 

Stecher, T.P. 1965, ApJ, 142, 1683 

Struve, O. 1929, MNRAS, 89, 567 

Struve, O., & Elver, C.T. 1938, ApJ, 88, 364 

Swings, R, & Rosenfeld, L. 1937, ApJ, 86, 483 

Tanaka, M., et al. 1996, PASJ, 48, L53 

Thompson, G.I., Nandy, K., Morgan, D.H., & Houziaux, L. 1988, MNRAS, 230, 429 
Toft, S., Hjorth, J., & Burud, 1. 2000, A&A, 357, 115 
Trumpler, R.J. 1930, PASP,42, 214 

Tuffs, R.J., Popescu, C.C., Velk, H.J., Kylafis, N.D., & Dopita, M.A. 2004, A&A, 419, 821 

Uchida, K.I., Sellgren, K., & Werner, M.W. 1998, ApJ, 493, L109 

Wang, J., Hall, P.B., Ge, J., Li, A., & Schneider, D.P. 2004, ApJ, 609, 589 

Weingartner, J.C., & Draine, B.T. 2001a, ApJ, 548, 296 

Weingartner, J.C., & Draine, B.T. 2001b, ApJS, 134, 263 

Whittet, D.C.B. 2003, Dust in the Galactic Environment (2nd ed; Bristol: loP) 

Witt, A.N., & Vijh, U.P. 2004, see WCD04, 115 

Witt, A.N., Bohlin, R.C., & Stecher, T.P. 1986, ApJ, 305, L23 

Witt, A.N., Clayton, G.C., & Draine, B.T, ed., 2004, ASP Conf. Ser. 309, Astrophysics of Dust 
(San Francisco: ASP) (hereafter WCD04) 

Witt, A.N., Gordon, K.D., & Furton, D.G. 1998, ApJ, 501, LI 11 

Witt, A.N., Smith, R.K., & Dwek, E. 2001, ApJ, 550, L201 

Witt, A.N., Lindell, R.S., Block, D.L., & Evans, R. 1994, ApJ, 427, 227 

Wolff, M.J., Clayton, G.C., Kim, S.H., Martin, P.G., & Anderson, C.M. 1997, ApJ, 478, 395 

Zubko, V.G., Dwek, E., & Arendt, R.G. 2004, ApJS, 152, 211 



